arXiv:1505.01497vl [astro-ph.GA] 6 May 2015 


Mon. Not. R. Astron. Soc. 000,[T|J^() Printed 8 May 2015 (MN IATeX style file v2.2) 

Comparative internal kinematics of the Hii regions in 
interacting and isolated galaxies: implications for massive 
star formation modes. 

Javier Zaragoza-Cardiel^’^*, John E. Beckman^Joan Font^’^, 

Begoha Garcia-Lorenzo^’^, Artemi Camps-Fariha^’^, Kambiz FathP, Philip A. James^, 
Santiago Erroz-Ferrer^’^, Jorge Barrera-Ballesteros^’^, Mauricio Cisternas^’^ 

^ Instituto de Astrofisica de Canarias, C/ Via Ldctea s/n, 38205 La Laguna, Tenerife, Spain 
‘^Department of Astrophysies, University of La Laguna, E-38200 La Laguna, Tenerife, Spain 
^CSIC, 28006 Madrid, Spain 

^Stoekholm Observatory, Department of Astronomy, Stoekholm University, AlbaNova, 106 91 Stockholm, Sweden 
^Astrophysics Research Institute, Liverpool John Moores University, I 46 Brownlow Hill, Liverpool L3 5RF, UK 


ABSTRACT 

We have observed 10 interacting galaxy pairs using the Fabry-Perot interferometer 
GHctFaS (Galaxy Hct Fabry-Perot system) on the 4.2m William Herschel Telescope 
(WHT) at the Observatorio del Roque de los Muchachos, La Palma. We present here 
the Hct surface brightness, velocity and velocity dispersion maps for the 10 systems 
we have not previously observed using this technique, as well as the physical prop¬ 
erties (sizes, Ha luminosities and velocity dispersion) of 1259 Hll regions from the 
full sample. We also derive the physical properties of 1054 Hll regions in a sample 
of 28 isolated galaxies observed with the same instrument in order to compare the 
two populations of Hll regions. We find a population of the brightest Hll regions for 
which the scaling relations, for example the relation between the Ha luminosity and 
the radius, are clearly distinct from the relations for the regions of lower luminosity. 
The regions in this bright population are more frequent in the interacting galaxies. 
We find that the turbulence, and also the star formation rate, are enhanced in the Hll 
regions in the interacting galaxies. We have also extracted the Ha equivalent widths 
for the Hll regions of both samples, and we have found that the distribution of Hll 
region ages coincides for the two samples of galaxies. We suggest that the SFR en¬ 
hancement is brought about by gas flows induced by the interactions, which give rise 
to gravitationally bound gas clouds which grow further by accretion from the flowing 
gas, producing conditions favourable to star formation. 

Key words: galaxies: interactions - stars: formation - galaxies: ISM - galaxies: 
kinematics and dynamics - (ISM:) H ii regions 


1 INTRODUCTION 


Galaxy mergers play an important role in galaxy evolution 
( White Rees|1978 ), but it is still not clear whether, and if 


Bournaud|2011 Tadhunter et al.|2011 

), whether they stim- 

ulate nuclear activity (jCanalizo et al. 

|2007| 

Bennert et al. 

20081 IGeorgakakis et al.|2009| Cisternas et a] 

.|2011| IRamos 


Almeida et al.|2011[ 2012| Bessiere et al. |2012| ), and what 

are the mechanisms by which they produce new structures 


E-mail:jzc@iac.es 


such as tails, bars and warps. They are important sites for 
the feedback processes, which are considered to be of con¬ 
siderable importance in producing realistic galaxies in the 


context of ACDM cosmological models. (Bournaud et al. 
2010 Hopkins et al.||2013 ) 


Because mergers imply distortions quantifiable most ac¬ 
curately by kinematical studies of the galaxies as a whole, 
3D studies of galaxy mergers are extremely important. Re¬ 


cently, Wisnioski et al. (2014) published a 3D study of galax¬ 


ies with redshift 0.7 < z < 2.7 claiming that the majority 
of galaxies are star forming and turbulent discs dominated 
by rotation and therefore not strongly affected by mergers. 
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However, nearby galaxy studies can give us better clues in 
practical terms about how to study star forming systems 
at higher redshifts, since we can achieve far better spatial 
resolution. 

Star formation enhancement in galaxy collisions is the¬ 
oretically well predicted for the central parts of the galaxies 
since collisions of galaxies tend to produce linear structures, 


interacting galaxies, Arp 270 and the Antennae galaxies. 


inducing strong gas inflows towards the nucleus (Mihos & 
Hernquist|1996 ). However, the off-nuclear peaks of star for¬ 


mation, such as are clearly present in the Antennae galaxies, 
are much less well understood, and they are likely to be at 
least as important ( Bournaud et al.|2010 Smith et al.|2014 ). 

Ka emission is one of the most useful indicators of mas¬ 
sive star formation in galaxies. Observing its full spectral 
profile is a very powerful tool for disentangling kinematics, 
but simply mapping the Ha surface brightness of a galaxy 
is already a comprehensive way to gain insight into global 


star formation. Kennicutt, Edgar, & Hodge (1989); Bradley 


et al. (2006) studied the properties of the Hii region popula¬ 


tions, and in particular their luminosities, hnding a break in 
the luminosity function at logLna = 38.6 dex, which some 


authors (Beckman et al. 2000) proposed could be used as 


a distance calibrator for galaxies. However, a physical basis 
for the break is desirable in order to strengthen its use as a 
distance calibrator, and while scenarios have been proposed 
( Beckman et al.|2000 Zaragoza-Cardiel et al.|2013 ), a clear 
explanation has not so far been agreed on. 


Terlevich & Melnick (1981) studied the scaling relations 


of extragalactic Hll regions measuring the sizes, luminosities, 
and velocity dispersions. They found that the luminosity, the 
size, and the velocity dispersion were correlated {L oc 
and L oc a^). The L — ay relation has also been proposed as 
a distance calibrator ( Terlevich Melnick||1981 Chavez et 


aL]|2014 ). For this it is desirable to agree on a physical ba¬ 
sis for the driver of the turbulence in Hll regions, but there 
are a number of possible processes, which may combine in 
different proportions to produce the observed line widths. 
The studies by ( Terlevich Melnick|1981 Zaragoza-Cardiel 


et al. 2013 20i4| Chavez et al. 2014) point towards self¬ 


gravity, at least in the more luminous regions, while 


Dyson 


(1979); Hippelein & Fried (1984); Rosa & Solf (1984); 


Chu 


& Kennicutt (1994) favour supernova explosions and stel¬ 


lar winds. Gutierrez et al. (2011) used HST data from M51 


to extract sizes and luminosities with unprecedented resolu¬ 
tion and found similar scaling relations to those cited in the 
Lila — R relation. 


However, our recent results (Zaragoza-Cardiel et al. 


2013) found different scaling relations for the brightest set 


of Hll regions in the interacting pair of galaxies Arp 270. We 
have found the same discrepancy in the Hii regions of the 


Antennae galaxies (Zaragoza-Cardiel et al. 2014) who ex¬ 


plain this as arising from two distinct populations of molecu¬ 
lar clouds as they find an equivalent dichotomy in the molec¬ 
ular cloud population of the Antennae. This change in the 
scaling relations changes the behaviour of the star formation 
rate (SFR, derived from Ha luminosity) dependency on the 
parameters (mass, radius) of the star forming regions. 

Here we present new observations of 10 interact¬ 
ing galaxies with the Fabry-Perot interferometer GHaFaS 
(Galaxy Ha Fabry-Perot System). We have measured the 
basic physical parameters (Lna, R, and ay) of each Hll re¬ 
gion and add equivalent data from the previous studies of 


Zaragoza-Cardiel et al. (2013 2014), respectively. We have 


extracted the same parameters for the Hll regions of 28 iso¬ 
lated galaxies from Erroz-Ferrer et al. (2015) in order to 


compare results from the two samples. 

For practical purposes we have called the galaxies which 
have no close neighbours and which therefore suffer negligi¬ 
ble interaction with other galaxies "isolated"even though the 
term is used in this article in a somewhat looser way than in 
studies dedicated to the careful dehnition of galaxies as iso¬ 
lated. The sample of isolated galaxies was selected to spread 
the whole morphological types, and different absolute mag¬ 
nitudes ranging from —20.6 to —17.7. 

In section §2 we present the galaxy sample, the obser¬ 
vations, and the data reduction. In section §3 we present 
the moment maps of the interacting galaxies. In section §4 
we explain how we derive the parameters of the Hil regions 
from those extracted directly from the observations. In sec¬ 
tion §5 we present the physical properties of the Hll regions 
in the interacting galaxies and in §6 those from the sample 
of isolated galaxies. In section §7 we discuss the differences 
between the two samples, making comparison with previous 
studies, and in section §8 we present a hypothesis to explain 
the zones of strong star formation well away from the galac¬ 
tic nuclei. Finally, in section §9 we present our conclusions. 


2 OBSERVATIONS AND GALAXY SAMPLE 


We observed the interacting galaxies with GHaFaS (Hernan- 
dez et ^|2008| ), a Fabry-Perot interferometer mounted on 
the WHT (William Herschel Telescope) at the Roque de los 
Muchachos Observatory, La Palma. The optical color com¬ 
posite images of the sample are in Fig. The observations 
were taken during several observing runs between December 
2010 and February 2013 (see Table [^. 

These galaxies were selected from the catalog of inter¬ 
acting galaxies by Arp (1966) and [Vorontsov-Velyamiii^ 
et al. ( 2003| (except NGC 2146 and UGC 3995, which are 
classihed as Galaxy pairs in de Vaucouleurs et al. (1976) 
and Karachentsev et al. ( 1985| ), respectively) according to 
the following criteria: 

(i) Nearby systems {Ysys < 10000 km/s), in order to have 
adequate spatial resolution. From ~ 30 pc to ~ 100 pc in our 
sample. 

(ii) Angular size which fits the GHaFaS FOV and allows 
satisfactory sky subtraction. 

(iii) Declination -20° < 6 < 80°, to ensure observations 
with air mass smaller than 1.6. 


GHaFaS has a circular FOV of 3.4 arcmin diameter. The 
etalon has a FSR (Free Spectral Range) of 8A in Ha, 
which corresponds to a 390km/s with a spectral resolution 
of 8km/s. The pixel size is 0.2 arcsec, which gave us the 
values of seeing limited angular resolution given in Table 
Further details about each galaxy observation are listed 
in Table The morphological type, apparent and absolute 
magnitude of the whole sample of interacting galaxies are in 
Table[^ The comparison of these properties with those of the 
sample of isolated galaxies from Erroz-Ferrer et al. (2015), 


yields that both of them spread over the whole morphologi¬ 
cal types, and that the interacting galaxies sample is brighter 
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Figure 1. Color composite images for the part of the sample of interacting galaxies we present here new observations. All of them are 
combined using SDSS data, except NGC 2146 (using DSS), and NGC 2992/3 (using data from du Pont telescope at Las Campanas 
observatory). 


in average, although less than an order of magnitude. Thus, 
the galaxy masses (derived from absolute magnitude) in the 
interacting galaxies are higher in average compared to the 
masses in the isolated galaxies. 


2.1 Data reduction 

GHoFaS is mounted at the Nasmyth focus of the WHT 
without an optical derotator. Correction for field rotation 


is therefore necessary, and was performed using the tech¬ 
nique described in detail in Blasco-Herrera et al. (2010). 
The observations are divided into individual cycles of 8 min, 
and each cycle is divided into 48 channels sepparated by 
8km/s and spanning the whole spectral range, which allows 
us to calibrate in velocity before applying digital derota¬ 
tion. We applied to the datacube for a given galaxy the 
procedures detailed in Daigle et al. (2006), including phase- 
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Name 

RA (J2000) 
(hh:mm:ss) 

Dec (J2000) 

(° ' ") 

Date 

dd/mm/yyyy 

Exp time 
(min) 

Seeing 

(arcsec) 

Dist 

(Mpc) 

NGC 520 

01:24:35.1 

+03:47:33 

18/11/2012 

181 

0.9 

34.4 

NGC 2146 

06:18:37.7 

+78:21:25 

23/12/2010 

181 

0.9 

21.9 

UGG 3995 

07:44:09.3 

+29:14:48 

19/11/2012 

181 

1.0 

60.6 

NGG 2782 

09:14:05.1 

+40:06:49 

17/11/2012 

181 

1.0 

23.4 

NGG 2992 

09:45:42.0 

-14:19:35 

19/11/2012 

156 

1.0 

31.6 

NGG 2993 

09:45:48.3 

-14:22:06 

18/11/2012 

148 

0.9 

30.5 

NGG 3769 

11:37:44.1 

+47:53:35 

29/05/2011 

156 

0.8 

16.1 

NGG 3786 

11:39:42.5 

+31:54:33 

28/02/2012 

205 

1.0 

50.9 

NGG 3788 

11:39:44.6 

+31:55:52 

28/02/2012 

205 

1.0 

50.1 

NGG 3991 

11:57:31.1 

+32:20:16 

04/02/2013 

181 

1.2 

43.7 


Table 1. Interacting galaxies observations log. 


Name 

rriB 

Mb 

morphology 

NGG 520 

12.2 

-20.48 

Sa 

NGG 2146 

10.53 

-21.17 

SBab 

UGG 3995 

13.52 

-20.39 

Sbc 

NGG 2782 

12.32 

-19.53 

SABa 

NGG 2992 

13.06 

-20.22 

Sa 

NGG 2993 

13.11 

-19.81 

Sa 

NGG 3395 (Arp 270) 

12.4 

-19.32 

Sc 

NGG 3396 (Arp 270) 

12.43 

-19.29 

SBm 

NGG 3769 

12.54 

-18.49 

Sb 

NGG 3786 

13.47 

-20.06 

SABa 

NGG 3788 

13.44 

-20.06 

SABa 

NGG 3991 

13.51 

-19.69 

Im 

NGG 4038 (Antennae) 

10.85 

-20.86 

SBm 

NGG 4039 (Antennae) 

11.04 

-20.67 

SBm 


Table 2. Whole interacting galaxies sample. The apparent mag¬ 
nitude and the morphology are taken from Hyperleda (http: 
//leda.univ-lyonl. fr |Paturel et ah] ( |2003| >). 


4 PHYSICAL PROPERTIES OF HII REGIONS 

We used astrodendrcQ a Python package to compute 
"dendrograms"of Astronomical data (Rosolowsky et al. 
2008) to extract the relevant parameters of the Hll regions 


from the datacubes. The parameters extracted are the Ha 
luminosity, Lna, the radius of the region, R, and the velocity 


dispersion, cf^ as explained in Rosolowsky & Leroy (2006); 


Rosolowsky et al. (2008), from where we have extracted a 


brief description in order to make it easier to understand the 
present article. 

We can describe a GHaFaS datacube as a collection of 
values, Ti, at a given point in the 3 dimensions which de¬ 
fine the coordinates of the datacube, two spatial dimension, 
[xi^yi)^ and the third the velocity along the line of sight, 
Vi. The method considers that a region is separated from 
the rest of the cube by an isosurface of value Tedge where 
T > Tedge inside the region and T < Tedge just outside the 
region. The method then estimates the major and minor 
axes as the mean values of the second spatial moments in 
the directions of these axes, respectively; 


correction, spectral smoothing, sky subtraction and gaussian 
spatial smoothing. 

The flux calibration of the GHaFaS datacubes is per¬ 
formed using a continuum-subtracted and flux-calibrated 
AG AM (Auxiliary-port GAMera, Benn et al. (2008)) Ha 
image. AG AM is an instrument mounted permanently at 
the WHT used for broad-band and narrow-band imaging. 
We have followed the procedure explained in |Erroz-Ferrer| 


et al. (2012), which consists in measuring the flux from a 


set of selected H H regions in both the GHaFaS cube and 
the AG AM image, and calibrating by direct comparison. 


3 MOMENT MAPS OF INTERACTING 
GALAXIES 


For each galaxy datacube, we have extracted the Ha surface 
brightness, velocity, and velocity dispersion maps using the 
tools described in Daigle et al. (2006). We show the maps in 
Fig. The complete sample of the moment maps is available 
through GDS. We have also included in the information sent 
to GDS the continuum subtracted Ha emission datacube. 


C^maj ,min 




EiTi 


( 1 ) 


where the sum is over the pixels inside the region, and 
X are the points lying on the major or the minor axis of the 
region. 

The equivalent radius is then, R = ?7^crmaj (Jmin where 
it is assumed that the region is spherical, so 77 = 1.91. 

The second moment is used to estimate the velocity 
dispersion along the velocity axis weighted by the datacube 
values 


(J V 


iEiV{vi-vr 


Er 


( 2 ) 


where 


y-region 


^^region ^ 


( 3 ) 


and the sum is over all the pixels inside the region de¬ 
fined by Tedge • Equation assumes that the observed Ha 


^ http://www.dendrograms.org 
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Figure 2. Moment maps of the interacting galaxies sample derived from the GHaFaS datacube. The colour versions of the maps are 
available in the electronic version of the article and the moment maps are available through CDS. 
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Figure 2 - continued 
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Figure 2 - continued 


profiles are gaussians. We plot a subsample of Ha profiles 
for a selection of the regions from our sample of galaxies 
observed with GHaFaS where we can check how Gaussian 
are the emission profiles in Fig. We might expect to 
find multiple components in the emission profiles as well 
as broad wings typical of Hll regions, however, the identifi¬ 
cation method differentiates between them if they are sep¬ 
arated in the datacube by a few times the rms, and selects 
the Gaussian core. A more detailed study of multiple com¬ 
ponents and expansive superbubbles using GHaFaS obser¬ 
vations is described in Gamps-Farina et al. (2015). 


The zeroth moment is used to estimate the flux F — 
J]]. Ti Sv SxSy, where Sx, Sy, and Sv are the pixel sizes. 

A key reason to use this method is that it is unbiased 
with respect to Tedge- We adopt Tedge = 4 x /rms, set the 
intervals used in the process of searching for different Hll 
regions to the same value as Tedge, and consider that the 
minimum area in pixels is the resolution of each observation. 
Three paradigms can be used with these method described 
in Rosolowsky et al. (2008). We use here the “bijection” 


paradigm. 

Although correcting Ha for dust attenuation would 
strengthen our conclusions, making this correction would 
not affect the scaling relation behaviour, as we showed for 
the Antennae galaxies ( Zaragoza-Gardiel et al.|2014 ) where 


there is quite heavy and variable dust attenuation (Brandi 
et al.|200'5 ). In a very detailed study using HST Ha and con¬ 


tinuum images of M51, Gutierrez & Beckman (2010) found 


that the dust attenuation is not a systematic function of Hll 
region luminosity, or size. This is explained using an inho¬ 


mogeneous model in which the bulk of the ionization occurs 
in denser clumps, in such a way that the dust attenuation 
depends on the mean clump size rather than on the overall 
size of the region. This implies that a statistical study of the 
type presented here should not be significantly affected. 

Given the Ha luminosity, Tna, and the size of the re¬ 
gion, R, we can derive the mean electron density, ne, fol¬ 


lowing Relano et al. (2005). Making a spherical first order 
approximation for the Hll regions ( Spitzer|1978 ): 


Lucx 
77 Rem 


= /^^^HaaH®a(i^O,T)2.46•10" 


■ nlRc 


(4) 


where hi/ua is the energy of an Ha photon, a^^{Ho,T) 
is the effective recombination coefficient of the Ha emission, 
and Rem is the radius in cm. Equation assumes no vari¬ 
ation of the filling factor and of the ionizing photon scape 
probability with the properties of the regions. However, a 
comprehensive study of the filling factor for the Hll regions 


in NGG 6946 by Gedres et al. (2013) implies that there 


should be little variation in the effect of the Filing factor 
on the mean electron density between larger and smaller 
Hll regions, while a detailed study by [Zurita et al. (2002) 
showed that the ionizing photon escape fraction does not 
show strong variations with luminosity (and size), as would 
be expected from regions with essentially clumpy structure 
( Giammanco et al.||2004 ) 

Equation [2 also assumes that the source of ionization is 
only ionizing photons emitted by the massive stars, and does 
not include other sources of ionization such as shocks. Previ- 
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Figure 3. Ho; profile for a sample of the brightest Hll regions (Region 1 in tables EUl and |A2p ) of a subsample of galaxies. 


ous results reported in Calzetti et al. (2004) for a selection of 
galaxies suggest that the proportion of ionization by shocks 
is rather small compared to the effect of photoionization. 
The diagnostics used to disentangle the two possible con¬ 
tributions to the ionization are in any case degenerate, but 
previous work gives an upper limiting value of 33% for the 
contribution due to the ionization by shocks, and a canonical 
value of 15% ( Hong et al.||2011 ). 

From Equation!^ we can derive the ionized gas mass, 
Mhii 


Mhii(M 0 ) = ^-KR^riemp = 1.57 x 10“^x (5) 

where H is in pc, Lua is in erg/s, and mp = 1.67 x 
10 “^^kg is the proton mass. 


As in Terlevich & Melnick ( 

19811; 

|Blasco-Herrera et al. 

(2010); Zaragoza-Gardiel et al. (2013 | 

|2014), the use of the 


velocity dispersion of the Hll region gives us really new infor¬ 
mation about the physics in those regions when compared to 
previous purely morphological studies ( Kennicutt, Edgar, 
Hodge|1989 Bradley et al.|20Q^ . We corrected the velocity 
dispersion for the instrumental, the natural = 3km/s), 
and the thermal (dy^th = 9.1km/s linewidths, assuming an 
isothermal Hil region temperature of T = lO^K), and sub¬ 
tracting them in quadrature from the observed width. The 
instrumental velocity dispersion is derived from the Neon 
calibration lamp using the emission line at 6598.9A. Eor each 
GHaEaS observation of one galaxy, a calibration datacube 
is taken, so we derive the instrumental velocity dispersion 



Figure 4. Emission profile of the Neon calibration lamp at 
6598.9A. 


fitting a Gaussian to the brightest pixel of the calibration 
datacube (see Eig. We only use as a valid velocity dis¬ 
persions those with line widths greater than 8 km/s, i.e., ve¬ 
locity dispersions greater than 4km/s, which is the effective 
velocity resolution of the instrument ( 8 km/s for GHaEaS). 

H 

We estimated the virial parameter, avir = 

( Bertoldi McK^|1992 ), using the corrected velocity dis¬ 
persion, the derived radius of the Hii regions, and using the 
mass of the ionized gas as the gas mass, Mgas = Mhii as 
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a first approximation. The stellar mass inside Hil regions is 
small (less than 4%) compared with the total Hil region mass 
( Relano et al.|200^ , so we use the size of the Hii region and 
not that of the stellar ionizing cluster to estimate the virial 
parameter. It is true that Hll regions, particularly large, lu¬ 
minous regions, have considerable neutral and molecular gas 
as well as ionized gas. In the larger more luminous regions, 
the ionized gas is an order of magnitude smaller than the 
neutral plus molecular gas ( Yang et al.|1996 Giammanco et 


al. 2004). The results presented in Zaragoza-Cardiel et al. 


(2014) suggest that a constant ionization fraction is a good 
approximation, at least, for the correlations we are going to 
study here. Thus, the real virial parameters for all the re¬ 
gions will be a nearly constant fraction of those estimated 
here, between 0.05 (Zaragoza-Cardiel et al. 2014) and 0.1 
(Yang et al. 1996 Giammanco et al. 2004), depending on 
the ionization fraction value. Since the virial parameter is in 
any case uncertain by a factor which depends in detail on 
the (non uniform) density distribution and deviations from 
sphericity ( [Bertoldi &: McKee|1992 ), the extra scatter in the 
uncertainty entailed by using ionized gas and a constant fac¬ 
tor to make up the total gas mass will not seriously affect 
our results. 


4.1 Uncertainties 


We have estimated the uncertainties in the measured pa¬ 
rameters following the bootstrapping method explained in 


Rosolowsky & Leroy (2006). The method estimates the un¬ 


certainties by randomly sampling the values of the datacube 
for each region allowing repeated values, and then subtracts 
the parameter (Luminosity, size or velocity dispersion) sev¬ 
eral times to estimate the standard deviation of these de¬ 
rived parameters. The uncertainty for each parameter is its 
standard deviation scaled up by the square root of the num¬ 
ber of pixels in one resolution element. This uncertainty does 
not include the intrinsic error of the flux in the data cubes. 


HII REGIONS IN THE INTERACTING 
GALAXIES 


We present the parameters of the Hii regions belonging to 
the interacting galaxies in Table [AT] 1 . We have included the 
data for the interacting pair Arp 270 presented in [Zaragoza-] 
Gardiel et al. (2013), where we have reanalyzed the Hll re¬ 


gion parameters using ASTRODENDRO package in substitu¬ 
tion of CLUMPFIND in order to produce a homogeneous study 
for the full set of galaxies. We present the reanalyzed Hll re¬ 
gions of Arp 270 in Table 1^. However, we have also used 
in the subsequent analysis the parameters of the Hll regions 
of the Antennae galaxies presented in jZaragoza-Gardiel etj 
al. (2014) in order to obtain a larger statistical sample of re¬ 


gions for the interacting galaxies. For these objects, we have 
1259 Hii regions in total. The parameter most sensitive to 
the noise level is the radius of the region, so in the follow¬ 
ing analysis we have removed the regions where the relative 
error in the radius is larger than 15%. After applying this 
exclusion level we were left with 537 Hll regions where the 
error in radius is sufficiently small, and for this sample the 
relative errors in the derived parameters: ionized gas mass, 
and electron density, are less than 20%. 


Name 

Al 

Li 

A 2 

L 2 

logRi 

NGC 2146 

2.4 

32.9 




Arp 244 

1.1 

35.3 

3.8 

30.2 

1.8 

NGC 520 

2.3 

32.6 

7.3 

21.5 

2.2 

UGC 3995 

2.9 

32.2 




NGC 3788 

1.8 

34.7 




NGC 3786 

3.1 

31.7 




NGC 2782 

1.3 

34.5 

3.6 

30.0 

1.9 

NGC 2993 

4.3 

29.0 




NGC 2992 

2.8 

32.5 




NGC 3991 

0.5 

35.9 

5.6 

25.0 

2.1 

Arp 270 

1.4 

34.8 

4.7 

28.4 

1.9 

NGC 3769 

1.4 

34.9 

5.4 

28.1 

1.7 


Table 3. Results of the single and double (when applicable) linear 
fits for interacting galaxies as defined in eq. 


5.1 Scaling relations of Hii regions 

In Fig. we plot the Lu^-R relations for the Hii regions for 
individual galaxies. We fit an x-axis error weighted single 
linear fit, and also a double continuous x-axis error weighted 
linear fit (4 free parameters) using non-linear least squares 
to fit the double linear fit function to the data. We weight 
the error on the x-axis since the radii of the regions are the 
strongest source of uncertainty. We have chosen between a 
continuous double linear fit or a single one depending on 
the Xred value, taking into account the addition of two free 
parameters more in the double linear fit. In the cases where 
Xred is smaller for the double linear fit than in the single 
one, we choose the double continuous x-axis error weighted 
linear fit: 


log(LHa) = Li + Al log(R); for log(R) < Ri 
log(LHa) = L 2 + A 2 log(R); for log(R) > Ri. 


(6) 


In the cases where x?ed i® smaller for the single linear 
fit than in the double one, we choose the single linear fit 
log(LHa) = Li + Al log(R). The results are in TableFor 
8 of the 12 galaxies in the sample, we found a regime where 
the exponent Ai and/or A 2 in the Lua — R relations is larger 
than three. We have found that for the galaxies with a dou¬ 
ble linear fit, the exponent A 2 is larger than three. Thus, 
the results for these larger (and brighter) Hll regions are 
different, when compared to the results of Terlevich Mel-| 
nick (1981); Gutierrez Beckman] (2010) for M51, where 


the maximum luminosity of its Hll regions is not as high 
as for our sample of interacting galaxies, and similar to re¬ 


cent results of Zaragoza-Gardiel et al. (2013 2014) for two 


of the systems included here, the Antennae and Arp 270. 
For two galaxies, NGG 2993, and NGG 3786 there are not 
enough fainter Hii regions to perform a double linear fit, al¬ 
though they show an exponent in the Lna — R relation also 
larger than three, but we can not separate them into two 
regimes. The Table and the Fig. are sorted by absolute 
magnitude, from brighter to fainter. The brightest regime 
where the exponent is larger than three is independent of 
the absolute magnitude. 

An exponent larger than 3 in the Lucx-R relation im¬ 
plies that the electron density (or the density of ionized gas) 
increases with the radius (and the luminosity) of the region. 
We show the Ue-R relation in Fig|^ (below each Lna — R 
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plot), where the electron density decreases with radius for 
the set of smaller (and fainter) Hil regions for the interact¬ 
ing galaxies. In those cases we find a two valued behaviour, 
density decreasing with radius in the small radius (i.e. low 
luminosity) part of the parameter space, and increasing with 
radius for large radii (i.e. large luminosities). This change is 
clearest for those systems with large numbers of luminous 
regions. 


5.2 Luoi-Uv envelope 


Since we have measurements of the velocity as a third di¬ 
mension, we have used this to deepen our study of the Hll re¬ 
gions. We plot the envelope in Fig|^(left) for all the 

Hii regions from the sample of interacting galaxies. |Relano| 
et al. (2005) suggested that the Hll regions on the envelope 
are virialized in the sense that they are gravitationally domi¬ 
nated since those regions are the ones showing the minimum 
velocity dispersion for a given luminosity, and the excess for 
the regions away from the envelope can be attributed to in¬ 
ternal motions which are not in quasi-equilibrium, such as 
recent contributions from stellar winds and supernovae. We 
have checked if the distance of an Hi I region, or a set of 
Hll regions, from the envelope could be used as a parameter 
to distinguish between two density regimes. We estimated 


the envelope as described in Relano et al. (2005), taking the 


Hll region with the minimum velocity dispersion for each 
luminosity bin, and fitting a linear relation. The result is 


log(Z/Haenv) — 2.1 X log((7^>) T 37.5 (7) 

where Lnaenv is in erg/s and ay is in km/s. We have 
divided the populations of Hil regions, those which are on 
the envelope and obey Lna > and those which are 

not on the envelope and do not obey the previous inequality. 
This condition is plotted as a dashed line in Fig[^(left). We 
have plotted the electron density versus radius for all the 
Hii regions from the sample of interacting galaxies in Fig|^ 
(right). The two populations are clearly separated in the 
rie-R relation plotted, where the grey points are for regions 
below the envelope, and the black points are on the envelope. 
The electron density decreases with the radius of the region 
for the regions below the envelope, and increases for the 
regions on the envelope, where the Ha luminosity depends 
super-linearly on the velocity dispersion (Equation]^. 


5.3 Virial parameter 


We have estimated the virial parameter, avir = 
which is the ratio between the kinetic and gravitational en¬ 
ergy ( Bertoldi McKee|1992 ), substituting the ionized gas 
mass, Mhii, for the total gas mass, as a practical first ap¬ 
proximation. The comparison between ionized and molec¬ 


ular gas in the Antennae galaxies (Zaragoza-Cardiel et al. 


2014) yields the result that the ionized gas mass must be 


a nearly constant fraction of the total gas mass. Since the 
virial parameter is affected by a further constant factor due 
to the density inhomogeneities and the non-sphericity of gas 
clouds, the use of the ionized gas mass yields a constant off¬ 
set for the virial parameter. As we are making a comparative 



log(FH«) (erg/s) 


(T^ R 

Figure 7. The virial parameter, Ovir = 5 versus Ha lumi- 

(_riV2gas 

nosity, I/Ha • We have used the mass of the ionized gas as the gas 
mass (in fact it should be a constant fraction of the total mass, 
but our qualitative conclusions do not depend on the value of this 
fraction, so we have not included it) Mgas = Afnii- 


analysis of the two sets of Hll regions, we will neglect the 
numerical values of these offsets in the first instance. 

We plot Gvir versus Lna in Fig|^ for all the Hii re¬ 
gions together. We have used a result from the literature 
which shows that the Hll regions from the "electron density 


decreasing "regime are pressure confined (Gutierrez et al. 


2011), and now we can check if the values of the virial pa¬ 
rameter agree with this, avir ^ 1 values, are for gas clumps 
confined by the external pressure. This explains why the 
gas density decreases with size since the external pressure 
decreases with the distance from the galactic plane, with a 
scale height comparable with the sizes of the Hll regions, so 
the bigger the region, the lower is the effective external pres¬ 
sure. Values of avir ~ 1 are for clouds where the self gravity 
is the dominant force. The fact that the more massive clouds 
take values of this order implies that for clouds with masses 
sufficiently large the bigger the region, the stronger is the 
gravitational field so the density is an increasing function 
of radius. In Fig|^we can see that for the Hll regions from 
the lower luminosity regime, the virial parameter lies in a 
range significantly greater than unity, as predicted for pres¬ 
sure confined Hll regions. In contrast for the Hll regions in 
the high luminosity regime, the virial parameter lies close 
to a nearly constant value with mass, whose value is, appar¬ 
ently different from unity. However, if we take into account 
the results from Zaragoza-Cardiel et al. (2014) which give a 
nearly constant ionized gas fraction of ~ 0.05, substituting 
the total gas mass for the ionized gas mass the virial param¬ 
eter does lie close to unity, implying that the brightest Hil 
regions are gravitationally bound. 


6 HII REGIONS IN ISOLATED GALAXIES 

We derived the parameters of the Hll regions in isolated 
galaxies using observational data from the sample of 28 iso¬ 
lated galaxies in Erroz-Ferrer et al. (2015) observed with 
the same instrument, GHaFaS, and on the same telescope, 
WHT, as our sample of interacting galaxies. We present the 
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log(i?)(pc) 


Figure 6. Left: Ha luminosity, versus velocity dispersion, (7^;, for Hll regions in interacting galaxies. We have plotted the fitted 

envelope as a solid line, and the displaced envelope condition explained in the paper, as a dashed line. Right: electron density, rie, versus 
the radius, R, for Hll regions on the envelope (black) and those regions below the envelope (grey). 


parameters in Table |A2p for 1054 Hll regions. Again, we 
have used in the following analysis of Hii regions in isolated 
galaxies, those regions with a fractional error in radius less 
than a 15% since the radius of the regions is the most sen¬ 
sitive parameter. After applying this criterion we retained 
1018 Hll regions from the isolated galaxy sample. 


6.1 Scaling relations 


We plot Lua versus R for the Hii regions in isolated galaxies 
in fig. 1^ 

Again, we have chosen between a continuous double lin¬ 
ear fit or a single one depending on the Xred value, taking 
into account the addition of two free parameters more in 
the double linear fit. The results of the linear fits are in Ta- 
bl^ We find an exponent Ni and/or N 2 larger than three 
in the Lua — R relation in only 6 of 28 isolated galaxies. 
Thus, the regime where the ionized gas density increases 
with size (and luminosity) in the isolated galaxies is much 
less frequent than in the interacting galaxies. 

In Fig. (below each Lua — R relation) we plot the 
electron density versus the radius of the Hll regions for each 
individual galaxy. For the regions in the isolated galaxies we 
find in some cases a regime where the electron density in¬ 
creases with radius but where this occurs there are relatively 
few Hii regions and the phenomenon is less frequent. Fig. [8 
and Tabl^ a re sorted by absolute magnitude (see |Erroz^ 
Ferrer et al. (2015)), from brighter to fainter. The regime 
where the exponent is larger than three is again indepen¬ 
dent on the absolute magnitude of the galaxy. 


Name 

Ni 

Li 

N2 

L2 

logRi 

NGC 3504 

2.6 

33.5 




NGC 5678 

2.5 

33.0 




NGG 2805 

3.1 

30.9 




NGC 5921 

2.2 

33.4 




NGC 6070 

3.1 

31.0 




NGC 4151 

1.7 

34.5 




NGC 864 

2.6 

32.5 




NGC 2543 

1.7 

34.5 




NGC 2748 

2.7 

32.2 




NGC 3041 

0.4 

36.6 

2.2 

33.0 

1.9 

NGC 2712 

1.7 

34.4 




NGC 5740 

2.4 

32.9 




NGC 6412 

2.6 

32.6 




NGC 6207 

0.8 

36.1 

3.0 

31.9 

1.8 

NGC 1073 

2.9 

32.5 




NGC 3423 

1.4 

34.3 

3.4 

30.6 

1.9 

NGC 428 

2.5 

33.3 




NGC 5334 

2.2 

33.8 




NGC 5112 

2.9 

32.2 




NGC 3403 

1.6 

34.4 




NGC 4639 

2.6 

32.2 




NGC 7241 

3.0 

31.8 




NGC 918 

2.4 

33.3 

7.4 

22.8 

2.0 

NGC 4324 

4.0 

30.2 




NGC 4389 

2.5 

32.5 




NGC 4498 

2.6 

32.8 




NGC 2541 

3.1 

31.5 

4.6 

28.6 

1.9 

NGC 2500 

2.6 

32.9 





Table 4. Results of the single and double (when applicable) linear 
fits for isolated galaxies, as defined in eq. 


6.2 Lucx-av envelope 

Following the same procedure as in the previous section in 
interacting galaxies, we plot the Ha luminosity versus ve¬ 
locity dispersion for Hll regions in isolated galaxies in Fig. 
I^(left). Again, we plot the envelope to the data 

( 8 ) 


implying that the SFR of the regions on the envelope 
depends super-linearly on the velocity dispersion. 

We plot the electron density versus the radius of the re¬ 
gions, separating those regions on the envelope (black) from 
those regions below the envelope (grey) in Fig[^ (right). For 
the set of isolated galaxies there is a deficiency of Hll re¬ 
gions with a big enough velocity dispersion to be measured 
reliably by our observations, compared with Hll regions in 


log(Z/Haenv) — 2.0 X log(cr^;) T 37.6 
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log(i?)(pc) log(i?)(pc) 



log(i?)(pc) log(i?)(pc) 


Figure 8 - continued 


interacting galaxies. This is clear when we compare Figs. 
(left) and(left). We conclude that the velocity dispersions 
in the Hll regions of interacting galaxies are larger presum¬ 
ably due to star formation feedback effects which enhance 
the dispersion. 

For the isolated galaxies we cannot detect the two 
regimes in the electron density-radius relation when we sep¬ 
arate the regions on and below the envelope because we do 
not have enough regions below the envelope (Figright), 
implying that the effects of star-formation feedback are sig¬ 
nificantly less than for the interacting galaxies, as we would 
expect. 


dominated by self-gravity rather than being pressure con¬ 
fined. However, we can see that there is clearly a smaller 
proportion of Hll regions in the brighter regime compared 
to the number of Hii regions in the brighter regime in inter¬ 
acting galaxies. 


7 COMPARISON 

We have studied the differences in the populations of Hii 
regions in the two samples of galaxies. We have plotted the 
Hii region luminosity functions for interacting galaxies and 
isolated galaxies, which are well represented by a power law 


6.3 Virial parameter 

We plot the virial parameter versus Ha luminosity in Fig. 

We see the same general behaviour as before, for the Hll 
regions in isolated galaxies the virial parameter decreases 
with luminosity, implying that the brightest Hll regions are 


N{L)dL = AL“dL 


(9) 


where N{L)dL is the number of Hll regions with lumi¬ 
nosity in the range [L, L + dL] . |Kennicutt, Edgar, Hodge 


(1989); Beckman et al.| ( |2000| ; [Bradley et al. (2006) showed 
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Figure 9. Left: Ha luminosity, versus velocity dispersion, ay, for Hll regions in isolated galaxies. We have plotted the fitted 

envelope as a solid line, and the displaced envelope condition explained in the text, as a dashed line. Right: Ue versus R for Hll regions 
on the envelope (black) and those regions below the envelope (grey). 



log(^HJ (erg/s) 


Figure 10. The virial parameter, avir = versus Ha 

Or iWgas 

luminosity, I/Ha- We have used for the gas mass, the ionized gas, 
Algas = AIhii- 


that for a galaxy with a sufficiently large population of Hll 
regions, the luminosity function is best fitted by a double 
power law. 

We plot the Hll region luminosity functions in Fig. mi 
showing in blue the Hil regions of the interacting galaxies, 
and in red those of the isolated galaxies, as well as the double 
power law fits. 

We have estimated the completeness in both samples 
dividing each sample into two subsamples characterized by 
the distance to the galaxies. Thus, we have two subsamples 
for two ranges of smaller and larger distances respectively, 
and we estimated the completeness limit where the slope of 
the luminosity function converge to the same value for both 
subsamples. The completeness limit for interacting galaxies 
is Lua — 37.3 dex while for isolated galaxies it is Lna = 
38 dex. 

For consistency with previous studies, we need to sub¬ 



Figure 11. Hll region luminosity function for interacting galax¬ 
ies (blue) and isolated galaxies (red). The completeness limit is 
drawn as a dashed vertical line in blue for interacting galaxies, 
and in red for isolated galaxies. The double fit is plotted as a 
solid line in blue for interacting galaxies, and in red for isolated 
galaxies. 

tract 1 from the slopes of the fits in Figj^to estimate the 
a value appearing in Equation since we have binned the 
data in logarithmic bins rather than linear bins. The results 
of the exponent in the luminosity function fit are 

a = (—1.17 ± 0.10); for log(LHa) < 38.6 dex 
a — (—1.71 zb 0.17); for log(LHa) > 38.6 dex 

for interacting galaxies, while for isolated galaxies the 
result is 

a = (—1.8 zb 0.2); for log(LHa) < 38.8 dex 
a = (-2.5 zb 0.3); for log(LHa) > 38.8 dex 

We draw the conclusion that the Hll region luminosity 
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functions for the two samples are different within the uncer¬ 
tainties. The Hil region population in the interacting galax¬ 
ies is generally brighter in Ha (and therefore the SFR(LHa) 
is higher) compared to isolated galaxies. The break in the 
luminosity function is compatible with those found in previ¬ 
ous studies ([Kennicutt, Edgar, k, Hodge|1989 Bradley et al. 


2006 IZaragoza-Cardiel et al.|2013 2014). Zaragoza-Cardiel 

et al.[(2013 2014) claimed that this double population of 


Hii regions is due to different scaling relations in the Lna-R 
relation. In this study the break in the scaling relation for 
different galaxies is not constant due to the scatter in the 
data from which the change in scaling relations is derived. 

We have applied the statistical Kolmogorov-Smirnov 
test to derive the probability that the Ha luminosity distri¬ 
butions we have obtained for interacting and isolated galax¬ 
ies are the same. The result is a probability of 10“^%. We 
can therefore conclude that the Ha luminosity distribution 
in the sample of interacting galaxies is different from that 
of the isolated galaxies sample, and as Fig. pT] shows, that 
the Hll regions interacting galaxies are brighter compared to 
those in isolated galaxies, because the luminosity function 
extends to higher values of the luminosity. 

We plot the distribution of rie in Fig(left). We can 
roughly compare the ionized gas density, rie, with the molec¬ 
ular gas density (Zaragoza-Cardiel et al. ]M4l l. The gas den¬ 
sity in the Hll regions from the interacting galaxies sample 
is clearly enhanced compared to those in the isolated galax¬ 
ies. Since the SFR depends super-linearly on the gas density, 
the enhancement in the gas density in the interacting galax¬ 
ies sample implies an strong enhancement in the SFR, and 
indeed of the star formation efficiency. 

We commented in the previous section that there is a 
deficiency of Hii regions with enough velocity dispersion to 
be measured by our observations. We also plot the distri¬ 
bution of (Jv in Fig^ (right) where we can see that the 
interaction of galaxies increases the turbulence in Hii re¬ 
gions. The two main mechanisms which can contribute to 
produce this effect are, firstly the continued accretion of gas 
onto the Hil regions even after the first massive stars have 
been formed, producing more massive regions, and there¬ 
fore increasing the turbulence. The other mechanism is the 
star formation feedback by stellar winds and supernova ex¬ 
plosions, implying that the feedback is stronger in galaxy 
interactions due to higher rates of star formation, and also 
the formation of more massive stars. 


7.1 Age distribution of Hll regions 

Using the Ha narrow band and continuum observations 
made with AC AM, we estimate the Ha equivalent widths, 
EW, of the Hii regions extracted in this study, in the same 
Zaragoza-Cardiel et al. ( 2013| 2014). Assuming 


way as m 

a direct relation between EW(Ha) and the age of the re¬ 
gion (Leitherer et al. 1999| ), we derive also the age of the 
Hii regions in this study, assuming solar metallicity and in- 
stantanous star formation. The results are included in 
Tables [AT] 1 and |Alp . We plot the normalized distribution of 
EWs and ages in Eig. (left and right respectively). The 
histograms for the interacting and the isolated galaxies coin¬ 
cide, so there is no differenciation of the age distributions of 
Hii regions in interacting and isolated galaxies. This could 
be taken as a point against our result of triggered star forma¬ 


tion in interacting galaxies since the Hii regions are forming 
on the same timescales for both samples of galaxies. How¬ 
ever, the age range in the distribution traced by Ha EW 
is much smaller than the typical merger time scale, ~ Gyr 
( Conseli^|2009 ). 


8 OFF-NUCLEAR PEAKS OF STAR 
FORMATION 

Galaxy interactions produce a reduction in axisymmetry. 


causing gas flows towards the central parts of galaxies (Keel 


et al.||1985 Mihos Hernquis^|1996 Bournaud||2011 ), giv¬ 

ing a well understood explanation of the nuclear starbursts 
observed in interacting galaxies. However, there are many 
examples of merging galaxies where there are intense star 
forming regions at sites away from the nuclei of the two 
galaxies ( [Smith et al.||20d7| Bournau'd||2011 ), such as NGG 
4038/9 (the Antennae galaxies), and NGG 4676 (the Mice). 

The two populations of Hll regions seem to be related to 
two underlying populations of GMGs which are converted to 
ionized gas having similar mass distributions as the placental 
molecular gas ( Zaragoza-Gardiel et al.|2QT4 ), at least in the 
more massive regions. 

We can imagine an inflow in a galaxy induced by the 
interaction with a companion galaxy, with velocity, Uin, to¬ 
wards the central parts, since galaxy interactions can pro¬ 
duce gas inflows. If a gas cloud is close to the inflow tra¬ 
jectory, the gas of the inflow will be accreted if the escape 
velocity, Uesc, from the gas cloud is larger than the velocity 
of the inflow, Vin. 

In order to estimate the escape velocity of the gas inflow 
from the gas cloud, we use a very simplistic model where we 
neglect the influence of the galaxy which is producing the 
inflow towards the center. A more realistic model in which 
the competing gravitational pulls of the cloud towards its 
centre and of the galaxy towards its centre would give more 
accurate results, but is outside the scope of the present ar¬ 
ticle. However, it is easy to show that for a major fraction 
of the original galactic disc the galactocentric force is sig¬ 
nificantly less than the perturbing force of individual clouds 
with radii greater than, say, 50 pc so the results presented 
here can be taken as qualitatively valid. The escape velocity 
from the surface of a gas cloud is 


2GM 


( 12 ) 


where b is the impact parameter of the inflow with re¬ 
spect to the centre of the cloud. Let us assume that the mass 
of a region is related to the radius by M(R) = Mq 


2GMoR^ N 

- - - = Vo R‘^ . 


(13) 


We have seen that for the Hll regions we can divide 
the scaling relations between luminosity and radius into 
two. The regions with high luminosity, significantly more 
abundant in the interacting galaxies, show exponents of 3 
or greater, while the regions with lower luminosity show 
exponents closer to 2. Previous work on molecular clouds 
has shown ( Roman-Duval et al.|2010 Lombardi et al.|201Q 
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Figure 12. Left: Distribution of electron densities Ue in Hll regions of interacting galaxies (blue) and isolated galaxies (red). Right: 
Distribution of velocity dispersion (jy in Hll regions of interacting galaxies (blue) and isolated galaxies (red). In both figures we plot two 
scales, one in blue at the right for interacting galaxies, and one at the left for isolated galaxies, in order to easily compare the different 
distributions. 
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Figure 13. Left: Normalized fistribution of Hll region equivalent widths, EW. Right: Normalized distribution of ages (under the 
assumption of instantanous star formation) derived using the equivalent width and Starburst99 model from |Leitherer et al.||T999| >. 


Kauffmann et al.|2010 Colombo et al.|2Q14 ) that molecular 


clouds in general have a relation between gas mass and cloud 
radius with an exponent close to 2, which is consistent with 
the value found by Larson (1981), frequently quoted as one 
of the Larson laws. For the interacting galaxy pair the An¬ 
tennae, Zaragoza-Cardiel et al. (2014) found an exponent of 


2.6 for the mass-radius relation of its most massive molecular 
clouds. It is straightforward to show that if the Larson law 
relating the velocity dispersion Gy to the cloud mass M©, 
{ay oc M°‘^) holds good, for a sequence of molecular clouds 
of increasing mass where the exponent in the mass-radius 
relation is 2 or greater, there will be a mass above which the 
clouds become gravitationally bound. In [Zaragoza-Cardiel] 


et al. (2014) we showed that this would account for the ten¬ 


dency of the virial parameter to approach unity for molecu¬ 
lar clouds of high enough mass (for the Antennae this implies 
cloud masses greater than some lO^.bM©). We also argued 
that any clouds with masses of this order in the galaxies 
prior to interaction would tend to accrete mass from the gas 
flows induced by the interaction. To pursue this argument 


further we can estimate, using somewhat simplified assump¬ 
tions, the escape velocity as a function of cloud size, taking 
two different scaling relations, one for the lower mass clouds, 
and the other for the high mass clouds. The corresponding 
relations can be expressed as: 


M = 3817i?^; for < 30pc . . 

M = 251i?^ ®^ fori? > 30pc 

We introduce Equation into Equation for differ¬ 
ent impact parameters obtaining the plots shown in Eig|14[ 
where we can clearly see the escape velocity enhancement 
in the case of the double mass-radius regime (solid lines) 
compared to the case of single mass-radius regime (dotted 
lines). 

An inflow of gas of Vin = 27km/s was observed for the 
interacting galaxy NGC 3396 ( jZaragoza-Cardiel et al. 12013 ). 
As an example, with this gas inflow velocity value we can 
conclude from figure that regions with radius larger than 
~ 60pc would accrete gas from the inflow in the double 
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Figure 14. Escape velocity, Vesc ,of the clouds as a function of 
their radius, for different values of the impact parameter, h. 
The solid lines represent the case of a double regime in the mass- 
radius relation (Equation |14| ) while the dotted lines represent the 
case of the single mass-radius relation with a canonical exponent 
oiN = 2. 


mass-radius regime, while only regions with radius larger 
than ~ 80pc would accrete gas in the single mass-radius 
regime. The latter may not be gravitationally bound prior 
to the merging process, but as they acquire more gas they 
become so. 

We suggest that our study, showing that the largest 
clouds are bound by their own gravitation, strongly supports 
the scenario where the gas accreted from the flowing gas 
passing near the these clouds, (which are already present in 
the isolated galaxies prior to interaction), is the cause of the 
star formation enhancement in interacting galaxies and can 
give rise to the off-nuclear peaks of star formation. Thus, 
although we still do not, at this stage, have a model for 
how the gravitationally dominated gas clouds form initially 
within the pre-merging galaxies, the inflow of gas in the 
interacting galaxies from the outer parts should increment 
the number of gravitationally dominated clouds. 


9 CONCLUSIONS 


We have presented new Fabry-Perot observations for 10 in¬ 
teracting galaxies: the Ha surface brightness, velocity dis¬ 
persion, and velocity dispersion maps. All of these values are 
available through CDS. We have amalgamated our tables of 
the properties of Hll regions for the sample of interacting 
galaxies presented here together with two interacting sys¬ 
tems, Arp 270, and the Antennae galaxies from our previous 


publications (Zaragoza-Cardiel et al. 


2013 


2014). We give 


here the catalogue of 1259 Hll regions, including the radius, 
Ha luminosity, velocity dispersion, electron density, ionized 
gas mass, the virial parameter, and the corresponding errors 
for each parameter. 

In order to compare the properties of Hll regions in in¬ 
teracting galaxies with those in isolated galaxies, we have 
included the parameters of Hii regions in a sample of 28 
isolated galaxies observed with the same instrument (|e rroz- 
Ferrer et al. 2015), obtaining these parameters (the same 
as those for the interacting galaxies) for 1054 Hii regions. 
We present a subsample of the brightest Hll regions of the 
catalogue in Tab 1 ^ A i'l l in for the interacting galaxies, and 
in Tabl ^A2^ for the isolated galaxies. The full catalogues 
of the Hii region samples are available as a machine read¬ 
able Tablein the electronic version of the article, as well as 
through CDS. 

The scaling relations, Lna — 7?, obtained here for in¬ 


teracting galaxies are different from previous studies for iso¬ 
lated galaxies ( Terlevich Melnick|198T Gutierrez & Beck- 


i| |2010 ) in most of the cases, and in agreement with re¬ 


cent results for the two interacting systems whose Hil regions 
we have included here, Zaragoza-Cardiel et al. (2013 2014). 
They differ in the exponent of the relation Lna cx: i? . We 
have found that for almost all the interacting galaxies, the 
Ana — R IS best represented by a double line, one for the 
smallest Hli regions (< lOOpc) similar to previous studies 
where the exponent N is smaller than three, and other for 
the largest Hll regions (> lOOpc) where the exponent N 
is larger than three, except for two galaxies where we do 
not have enough fainter Hll regions to fit a double linear 
fit and we find that the exponent N is larger than three, 
independently of the radius. If >3, then the electron (or 
the ionized gas mass) density derived from equations and 
increases with size, in contrast with previous results in 
isolated galaxies where the electron density in Hii regions 
decreases with size ( Terlevich Melnick||1981 Gutierrez & 
Beckma'n||2010 ). 


In the sample of isolated galaxies, we have found that 
the scaling relation Lna — R is similar to that found in pre¬ 
vious studies in 22 of our 28 isolated galaxy sample, with 
slopes smaller than three. 

The differences in the exponents of the Lhq; — R relation 
are independent of the absolute magnitude of the galaxies. 

The Lua — cFv plot shows an envelope which has been de¬ 
tected in previous published work, whose authors suggested 


that the regions on this envelope are virialized (Relano et al. 


2005). In the case of interacting galaxies, the regions on the 


envelope are the ones where the gas density increases with 
size. In the case of isolated galaxies we do not have suffi¬ 
cient regions with large enough velocity dispersion to reach 
a clear conclusion. The Ha luminosity, and the derived star 
formation rate of the Hll regions on the envelope depend 
super-linearly on the velocity dispersion. We have obtained 
the virial parameter avir which is an estimate of the ratio 


between the kinetic and the gravitational energy (Bertoldi & 
McKee||1992 ), and have used this to show that the brighter 


the region, the more clearly it is held together by self gravity 
rather than external pressure. 

We have quantihed the differences between the two pop¬ 
ulations of Hii regions, those from interacting galaxies, and 
those from the isolated galaxies. Comparing the luminos¬ 
ity functions we hnd that the Hll regions in the interacting 
galaxies are on average brighter than the Hll regions from 
the isolated galaxies. The K-S test performed on the two dis¬ 
tributions of Ha luminosities gives a probability of 10“^% 
for the hypothesis that the two populations of Hll regions 
are really the same. 

The histograms of electron density, and velocity disper¬ 
sion show that on average the Hll regions from interacting 
galaxies are denser and more turbulent than those from iso¬ 
lated galaxies. Since the star formation rate depends super- 
linearly on gas density, the enhancement in the gas density 
in interacting galaxies implies a signihcant enhancement in 
the star formation efficiency. 

However the distributions of equivalent widths and ages 
of Hii region coincide for the two samples of galaxies, show¬ 
ing that the formation of the young massive stellar compo¬ 
nent traced by Ha extends over a similar range in time for 
interacting and isolated galaxies. 
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However the sample of interacting galaxies is, on aver¬ 
age, brighter than that of the isolated galaxies. We could 
conclude that the presence of the population of brighter Hll 
regions in the interacting galaxies might be simply due to 
their greater masses. Nevertheless, the number of Hll regions 
in an interacting galaxy is bigger than the number in an iso¬ 
lated galaxy with the same absolute magnitude, suggesting 
that interactions do in fact increase the star formation rate. 

Based on the results from Wei et al. (2012) and 


Zaragoza-Cardiel et al. (2014), we claim that the two pop¬ 


ulations of Hil regions are related to two populations of 
GMC’s. Romeo & Agertz (2014) claim that Toomre stabil¬ 
ity criteria depends on the values of exponents in the Larson 
laws. In the case of the exponent in the mass-size relation, 
largest values of this exponent implies that larger scales are 
unstable, in agreement with the presence of larger star form¬ 
ing regions in the high mass regime. The presence of massive 
clumps of star formation in nearby galaxies as those studied 
here can improve our understanding about the clumpy star 


forming discs observed at higher redshifts (Elmegreen et al. 


2007). 


We claim that the brightest and most massive star form¬ 
ing regions can accrete gas that is inflowing towards the cen¬ 
tral parts of the galaxies, or has suffered a perturbation in its 
original orbit, induced by interactions. |Blasco-Herrera et ah] 
(2013) found in their sample of 11 starburst galaxies that 8 
of them show evidence of a recent merger while they are still 
dominated by rotation, in agreement with a picture where 
the interactions triggers flows and star formation without 
messing up the host galaxy’s kinematics. This could be a 
plausible scenario for the off-nuclear peaks of star forma¬ 
tion produced in galaxy collisions, a scenario which requires 
further observational and theoretical exploration. 


ACKNOWLEDGMENTS 


Based on observations made with the William Herschel Tele¬ 
scope operated on the island of La Palma by the Isaac 
Newton Group of Telescopes in the Spanish Observato- 
rio del Roque de los Muchachos of the Institute de As- 
troffsica de Ganarias. This research has been supported 
by the Spanish Ministry of Economy and Gompetitive- 
ness (MINEGO) under the grants AYA2007-67625-GO2- 
01, AYA2009-12903 and AYA2012-39408-G02-02. JEB ac¬ 
knowledge financial support to the DAGAL network from 
the People Programme (Marie Gurie Actions) of the Eu¬ 
ropean Union’s Seventh Eramework Programme EP7/2007- 
2013/ under REA grant agreement number PITN-GA-2011- 
289313. 


This research made use of Astropy, a community- 


developed core Python package for Astronomy ( Astropy 
Gollaboration et al.|2013 ), APLpy, an open-source plotting 


package for Python hosted at http://aplpy.github.com, and 
ASTRODENDRO, a Python package to compute dendrograms 
of Astronomical data (http://www.dendrograms.org/). We 
acknowledge the usage of the HyperLeda database 
(http: //leda.univ-lyonl .fr). 

We thank Roberto Terlevich, the referee of the paper, 
whose comments have led to important improvements on the 
original version of the paper. 


REFERENCES 

Arp, H. 1966, ApJS, 14, 1 

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et 
al. 2013, A&A, 558, A33 

Beckman, J. E., Rozas, M., Zurita, A., Watson, R. A., & 
Knapen, J. H. 2000, AJ, 119, 2728 
Benn, C., Dee, K., Sz Agocs, T. 2008, Proc. SPIE, 7014, 
70146X 

Bennert, N., Canalizo, G., Jungwiert, B., Stockton, A., 
Schweizer, E., Peng, C. Y., Lacy, M. 2008, ApJ, 677, 846 
Bessiere, P. S., Tadhunter, C. N., Ramos Almeida, C., 
Villar-Martm, M. 2012, MNRAS, 426, 276 
Bertoldi, E., & McKee, C. E. 1992, ApJ, 395, 140 
Blasco-Herrera, J., Eathi, K., Beckman, J., et al. 2010, MN¬ 
RAS, 407, 2519 

Blasco-Herrera, J., Eathi, K., Ostlin, G., Eont, J., & Beck¬ 
man, J. E. 2013, MNRAS, 435, 1958 
Bournaud, E., Elmegreen, B. G., Teyssier, R., Block, D. L., 
& Puerari, 1. 2010, MNRAS, 409, 1088 
Bournaud E., 2011, EAS, 51, 107 

Bradley T. R., Knapen J. H., Beckman J. E., Eolkes S. L., 
2006, A&A, 459, L13 

Brandi, B. R., Clark, D. M., Eikenberry, S. S., et al. 2005, 
ApJ, 635, 280 

Calzetti, D., Harris, J., Gallagher, J. S., HI, et al. 2004, 
AJ, 127, 1405 

Camps-Earina, A., Zaragoza-Cardiel, J., Beckman, J. E., 
et al. 2015, MNRAS, 447, 3840 
Canalizo, G., Bennert, N., Jungwiert, B., Stockton, A., 
Schweizer, E., Lacy, M., Peng, C. 2007, ApJ, 669, 801 
Cedres, B., Beckman, J. E., Bongiovanni, A., et al. 2013, 
ApJL, 765, LL24 

Chavez, R., Terlevich, R., Terlevich, E., et al. 2014, MN¬ 
RAS, 442, 3565 

Chu, Y.-H., & Kennicutt, R. C., Jr. 1994, ApJ, 425, 720 
Cisternas M., et ah, 2011, ApJ, 726, 57 
Colombo, D., Hughes, A., Schinnerer, E., et al. 2014, ApJ, 
784, 3 

Conselice, C. J. 2009, MNRAS, 399, L16 
Daigle, O., Carignan, C., Hernandez, O., Chemin, L., & 
Amram, P. 2006, MNRAS, 368, 1016 
Dyson, J. E. 1979, A&A, 73, 132 

Elmegreen, D. M., Elmegreen, B. G., Ravindranath, S., & 
Coe, D. A. 2007, ApJ, 658, 763 
Erroz-Eerrer, S., Knapen, J. H., Eont, J., et al. 2012, MN¬ 
RAS, 427, 2938 

Erroz-Eerrer, S., Knapen, J. H., Leaman, R., et al. 2015, 
arXiv:1504.06282 

Georgakakis, A., et al. 2009, MNRAS, 397, 623 
Giammanco, C., Beckman, J. E., Zurita, A., & Relano, M. 
2004, A&A, 424, 877 

Gutierrez, L., & Beckman, J. E. 2010, ApJL, 710, L44 
Gutierrez, L., & Beckman, J. 2010, The Impact of HST on 
European Astronomy, 77 

Gutierrez, L., Beckman, J. E., & Buenrostro, V. 2011, AJ, 
141, 113 

Hernandez, O., Eathi, K., Carignan, C., et al. 2008, PASP, 
120, 665 

Hippelein, H., Sz Pried, J. W. 1984, A&A, 141, 49 
Hong, S., Calzetti, D., Dopita, M. A., et al. 2011, ApJ, 731, 
45 


















24 Zaragoza-Cardiel et al. 


Hopkins, P. F., Keres, D., Onorbe, J., et al. 2013, 
arXiv:1311.2073 

Karachentsev, L, Lebedev, V., & Shcherbanovski, A. 1985, 
Bulletin d’Information du Centre de Donnees Stellaires, 
29, 87 

Kauffmann, J., Pillai, T., Shetty, R., Myers, P. C., & Good¬ 
man, A. A. 2010, ApJ, 716, 433 
Keel, W. C., Kennicutt, R. C., Jr., Hummel, E., & van der 
Hulst, J. M. 1985, AJ, 90, 708 
Kennicutt R. C., Jr., Edgar B. K., Hodge P. W., 1989, ApJ, 
337, 761 

Larson R. B., 1981, MNRAS, 194, 809 
Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, 
ApJS, 123, 3 

Lombardi, M., Alves, J., & Lada, C. J. 2010, A&A, 519, 
L7 

Mihos, J. C., & Hernquist, L. 1996, ApJ, 464, 641 
Moiseev, A. V., Tikhonov, A. V., & Klypin, A. 2014, 
arXiv:1405.5731 

Paturel G., Petit C., Prugniel P., Theureau G., Rousseau 
J., Brouty M., Dubois P., Cambresy L., 2003, A&A, 412, 
45 

Ramos Almeida, C., Tadhunter, C. N., Inskip, K. J., Mor- 
ganti, R., Holt, J., Dicken, D. 2011a, MNRAS, 410, 1550 
Ramos Almeida, C., et al. 2012, MNRAS, 419, 687 
Relano, M., Beckman, J. E., Zurita, A., Rozas, M., & Gi- 
ammanco, C. 2005, A&A, 431, 235 
Roman-Duval, J., Jackson, J. M., Heyer, M., Rathborne, 
J., & Simon, R. 2010, ApJ, 723, 492 
Romeo, A. B., k Agertz, O. 2014, MNRAS, 442, 1230 
Rosa, M., k Solf, J. 1984, A&A, 130, 29 
Rosolowsky, E., k Leroy, A. 2006, PASP, 118, 590 
Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., k Good¬ 
man, A. A. 2008, ApJ, 679, 1338 

Smith, B. J., Struck, C., Hancock, M., et al. 2007, AJ, 133, 
791 

Smith, B. J., Soria, R., Struck, C., et al. 2014, AJ, 147, 60 
Somerville, R. S., Hopkins, P. E., Cox, T. J., Robertson, 
B. E., Hernquist, L. 2008, MNRAS, 391, 481 
Spitzer L., 1978, ppim.book, 

Tadhunter, C., et al. MNRAS, 412, 960 
Terlevich, R., k Melnick, J. 1981, MNRAS, 195, 839 
Toomre, A. 1977, Evolution of Galaxies and Stellar Popu¬ 
lations, 401 

de Vaucouleurs, G., de Vaucouleurs, A., k Corwin, J. R. 
1976, Second reference catalogue of bright galaxies, 1976, 
Austin: University of Texas Press., 0 
Vorontsov-Velyaminov, B. A., Noskova, R. L, k Arkhipova, 
V. P. 2003, VizieR Online Data Catalog, 7236, 0 
Wei, L. H., Keto, E., k Ho, L. C. 2012, ApJ, 750, 136 
White, S. D. M., k Rees, M. J. 1978, MNRAS, 183, 341 
Wisnioski, E., Eorster Schreiber, N. M., Wuyts, S., et al. 
2014, arXiv:1409.6791 

Yang, H., Chu, Y.-H., Skillman, E. D., k Terlevich, R. 
1996, AJ, 112, 146 

Zaragoza-Cardiel, J., Eont-Serra, J., Beckman, J. E., et al. 
2013, MNRAS, 432, 998 

Zaragoza-Cardiel, J., Eont, J., Beckman, J. E., et al. 2014, 
MNRAS, 445, 1412 

Zurita, A., Beckman, J. E., Rozas, M., k Ryder, S. 2002, 
AVA, 386, 801 


APPENDIX A: PROPERTIES OF THE HII 
REGIONS 



Internal kinematics of the Hll regions in interacting and isolated galaxies 25 


Table Al. Physical properties of the brightest Hll regions in interacting galaxies derived as described in section §4. The whole Table is 
available as a machine readable Tablein the electronic version of the paper and through CDS. 


N 

Galaxy 

RA 

Dec 

log(LHa) 

Rmi 

log(MHii) 

PHII 

(7 X) 

e^vir 

EWh„(A) 

log(Age)(yr 



hh:mm:ss 

° / // 

erg/s 

pc 

M© 

M0/pc3 

km/s 




1 

NGC2146 

6:18:22.7 

78:22:25.6 

39.04 ±0.04 

114 ±4 

5.9 ±0.06 

0.127 ±0.008 

15.0 ± 1.0 

40.0 ±9.0 

39.9 

6.82 

2 

NGG2146 

6:18:31.1 

78:21:39.6 

38.83 ±0.04 

116 ±4 

5.8 ±0.07 

0.097 ±0.007 

17.0 ± 1.0 

60.0 ± 10.0 

117.71 

6.78 

1 

Arp244 

12:1:55.6 

-18:52:52.0 

40.2 ± 0.03 

280 ± 5 

7.06 ±0.05 

0.126 ±0.006 

32.8 ±0.7 

30.0 ±3.0 

202.62 

6.76 

2 

Arp244 

12:1:53.9 

-18:52:42.7 

39.97 ±0.03 

253 ±6 

6.88 ±0.05 

0.111 ±0.006 

19.1 ±0.7 

14.0 ±2.0 

63.77 

6.8 

1 

NGG520 

1:24:35.6 

3:47:19.0 

39.12 ±0.04 

255 ±9 

6.46 ±0.07 

0.042 ± 0.003 

14.1 ±0.8 

20.0 ±4.0 

2.1 

7.19 

2 

NGG520 

1:24:35.5 

3:47:30.2 

38.57 ±0.03 

187 ±8 

5.99 ±0.08 

0.035 ± 0.003 

< 4.0 


10.93 

7.05 

1 

UGG3995 

7:44:11.3 

29:15:13.1 

40.29 ±0.04 

270 ± 10 

7.09 ±0.08 

0.14 ±0.01 

43.0 ±3.0 

50.0 ± 10.0 

17.66 

6.95 

2 

UGG3995 

7:44:11.1 

29:14:34.6 

40.06 ± 0.06 

390 ± 30 

7.2 ±0.1 

0.064 ± 0.009 

8.0 ± 1.0 

1.7±0.9 

11.48 

7.05 

1 

NGG3788 

11:39:42.5 

31:54:34.0 

39.7 ±0.02 

280 ± 10 

6.81 ±0.07 

0.072 ± 0.005 

11.9 ±0.9 

7.0 ±2.0 

23.29 

6.88 

2 

NGG3788 

11:39:42.3 

31:54:32.6 

39.58 ±0.04 

200 ± 20 

6.5 ±0.1 

0.1 ±0.01 

< 4.0 


14.92 

7.02 

1 

NGG3786 

11:39:45.0 

31:55:30.7 

39.86 ±0.02 

300 ± 10 

6.93 ±0.08 

0.078 ± 0.006 

16.5 ±0.7 

11.0 ±2.0 

18.88 

6.88 

2 

NGG3786 

11:39:45.3 

31:55:24.0 

39.69 ±0.02 

240 ± 10 

6.71 ±0.09 

0.088 ± 0.008 

7.0 ± 1.0 

3.0 ± 1.0 

8.79 

7.09 

1 

NGG2782 

9:14:5.1 

40:6:43.8 

39.8 ±0.02 

166 ±6 

6.52 ±0.07 

0.17 ±0.01 

28.0 ±3.0 

40.0 ± 10.0 

63.41 

6.8 

2 

NGG2782 

9:14:4.1 

40:6:47.2 

39.34 ±0.01 

128 ±5 

6.12 ±0.06 

0.15 ±0.01 

21.0 ±2.0 

50.0 ± 10.0 

21.63 

6.96 

1 

NGG2993 

9:45:48.2 

-14:22:7.8 

40.07 ±0.01 

258 ± 7 

6.94 ±0.04 

0.121 ±0.005 

12.0 ± 1.0 

5.0 ± 1.0 

135.88 

6.78 

2 

NGG2993 

9:45:50.8 

-14:23:48.2 

39.21 ±0.02 

175 ±9 

6.26 ±0.09 

0.081 ± 0.007 

< 4.0 


0.92 

7.31 

1 

NGG2992 

9:45:42.5 

-14:19:30.5 

39.18 ±0.05 

100 ±8 

5.9 ±0.1 

0.18 ±0.03 

16.0 ±3.0 

40.0 ±20.0 

10.78 

7.05 

2 

NGG2992 

9:45:42.2 

-14:19:22.6 

39.16 ±0.04 

127 ±9 

6.0 ±0.1 

0.12 ±0.02 

6.0 ±2.0 

5.0 ±4.0 

18.63 

6.96 

1 

NGG3991 

11:57:29.6 

32:19:49.4 

40.44 ±0.02 

470 ±8 

7.52 ±0.03 

0.076 ± 0.003 

20.9 ±0.8 

7.2 ±0.9 

47.72 

6.81 

2 

NGG3991 

11:57:29.6 

32:19:51.7 

39.47 ±0.02 

350 ± 10 

6.84 ±0.06 

0.039 ± 0.002 

8.0 ±0.7 

< 4.0 ± 1.0 

68.51 

6.8 

1 

Arp270 

10:49:54.1 

32:59:44.4 

40.22 ±0.04 

229 ±4 

6.94 ±0.05 

0.175 ±0.009 

20.4 ±0.9 

13.0 ±2.0 

2.43 

7.17 

2 

Arp270 

10:49:50.4 

32:59:2.4 

40.15 ±0.04 

146 ± 2 

6.62 ±0.05 

0.31 ±0.01 

35.0 ± 1.0 

49.0 ± 7.0 

96.73 

6.79 

1 

NGG3769 

11:37:51.8 

47:52:46.4 

38.31 ±0.05 

66 ±3 

5.17 ±0.09 

0.12 ±0.01 

< 4.0 


227.02 

6.76 

2 

NGG3769 

11:37:51.0 

47:52:51.5 

38.14 ±0.05 

70 ±4 

5.1 ±0.1 

0.09 ±0.01 

< 4.0 


148.42 

6.77 
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Table A2. Physical properties of the brightest Hll regions in isolated galaxies derived as described in section §4. The whole Table is 
available as a machine readable Tablein the electronic version of the paper and through CDS. 


N 

Galaxy 

RA 

Dec 

log(AHa) 

^HII 

log(MHii) 

PHII 

Cv 

ttvir 

EWHa(A) 

log(Age)(yr) 



hh:mm:ss 

° / // 

erg/s 

pc 

M© 

Mq/pc^ 

km/s 




1 

NGC3504 

11:3:10.4 

27:58:59.4 

39.61 ± 0.04 

162 ±5 

6.41 

± 

0.06 

0.145 ±0.009 

14.4 ±0.9 

15.0 ±3.0 

17.03 

6.98 

2 

NGG3504 

11:3:10.0 

27:58:47.2 

39.3 ±0.01 

80 ±7 

5.8 

± 

0.1 

0.29 ±0.04 

6.0 ±3.0 

6.0 ±7.0 

17.89 

6.96 

1 

NGG5678 

14:32:12.6 

57:53:49.0 

40.07 ±0.02 

430 ±9 

7.28 

± 

0.04 

0.057 ±0.002 

13.3 ±0.3 

4.7 ±0.5 

29.49 

6.83 

2 

NGG5678 

14:32:5.8 

57:55:6.5 

40.0 ±0.03 

380 ± 10 

7.15 

± 

0.06 

0.064 ± 0.004 

13.0 ± 1.0 

6.0 ± 1.0 

5.58 

7.1 

1 

NGG5921 

15:21:56.6 

5:4:12.4 

39.22 ±0.02 

320 ± 10 

6.67 

± 

0.06 

0.033 ± 0.002 

< 4.0 


1.99 

7.2 

2 

NGG5921 

15:21:53.6 

5:3:22.3 

39.15 ±0.03 

289 ±9 

6.56 

± 

0.06 

0.036 ± 0.002 

< 4.0 


16.31 

7.0 

1 

NGG6070 

16:9:55.2 

0:42:3.3 

39.57 ±0.04 

331 ±9 

6.86 

± 

0.06 

0.047 ±0.003 

18.0 ±0.9 

17.0 ±3.0 

79.3 

6.8 

2 

NGG6070 

16:9:56.4 

0:42:27.8 

39.45 ± 0.04 

279 ±9 

6.68 

± 

0.07 

0.053 ± 0.004 

10.4 ±0.7 

7.0 ±2.0 

120.82 

6.78 

1 

NGG4151 

12:10:32.6 

39:24:21.1 

38.21 ±0.01 

148 ±5 

5.65 

± 

0.06 

0.033 ± 0.002 

< 4.0 


126.03 

6.78 

2 

NGG4151 

12:10:32.6 

39:24:20.8 

38.18 ±0.02 

107 ±6 

5.43 

± 

0.09 

0.052 ± 0.005 

< 4.0 


138.64 

6.78 

1 

NGG864 

2:15:27.6 

6:0:8.8 

40.01 ±0.03 

204 ±4 

6.76 

± 

0.04 

0.162 ±0.007 

27.2 ±0.8 

30.0 ±4.0 

54.09 

6.81 

2 

NGG864 

2:15:28.7 

6:0:43.3 

39.78 ± 0.04 

291 ±6 

6.88 

± 

0.05 

0.073 ± 0.004 

17.5 ±0.7 

14.0 ±2.0 

27.61 

6.83 

1 

NGG2543 

8:12:56.2 

36:14:42.7 

39.13 ±0.04 

220 ± 20 

6.4 

± 

0.1 

0.052 ± 0.007 

< 4.0 


36.42 

6.82 

2 

NGG2543 

8:12:57.2 

36:14:48.6 

39.07 ±0.02 

320 ± 10 

6.59 

± 

0.08 

0.027 ±0.002 

< 4.0 


32.49 

6.82 

1 

NGG2748 

9:13:32.6 

76:27:55.0 

39.33 ± 0.03 

292 ±7 

6.66 

± 

0.05 

0.043 ± 0.002 

11.0 ±0.8 

9.0 ±2.0 

9.14 

7.09 

2 

NGG2748 

9:13:30.6 

76:27:41.9 

39.23 ± 0.03 

270 ±9 

6.55 

± 

0.07 

0.043 ± 0.003 

6.0 ±0.7 

3.0 ± 1.0 

9.95 

7.07 

1 

NGG3041 

9:53:2.4 

16:39:41.7 

39.02 ± 0.03 

293 ±7 

6.5 ±0.06 

0.03 ± 0.002 

4.5 ±0.7 

2.2 ±0.8 

22.31 

6.93 

2 

NGG3041 

9:53:6.5 

16:41:35.4 

39.02 ± 0.02 

381 ±7 

6.67 

± 

0.03 

0.0202 ± 0.0007 

< 4.0 


1.4 

7.25 

1 

NGG2712 

8:59:33.4 

44:55:27.5 

39.65 ± 0.04 

230 ± 10 

6.7 

± 

0.1 

0.09 ± 0.009 

12.2 ±0.8 

9.0 ±2.0 

51.86 

6.81 

2 

NGG2712 

8:59:32.7 

44:55:14.4 

39.21 ±0.03 

300 ±8 

6.61 

± 

0.05 

0.036 ± 0.002 

11.9 ±0.8 

12.0 ±3.0 

38.87 

6.82 

1 

NGG5740 

14:44:23.7 

1:40:52.4 

39.33 ± 0.02 

185 ±5 

6.35 

± 

0.05 

0.086 ± 0.004 

42.0 ±2.0 

170.0 ± 30.0 

39.41 

6.82 

2 

NGG5740 

14:44:25.3 

1:40:51.6 

39.18 ±0.02 

361 ±9 

6.72 

± 

0.04 

0.026 ±0.001 

5.6 ±0.5 

2.5 ±0.6 

15.91 

6.92 

1 

NGG6412 

17:29:27.3 

75:41:25.2 

39.96 ± 0.03 

404 ±7 

7.18 

± 

0.04 

0.055 ± 0.002 

19.3 ±0.7 

12.0 ±2.0 

2.48 

7.17 

2 

NGG6412 

17:29:27.4 

75:41:25.4 

39.33 ± 0.02 

340 ± 10 

6.75 

± 

0.06 

0.035 ± 0.002 

6.4 ±0.9 

3.0 ± 1.0 

2.39 

7.17 

1 

NGG6207 

16:43:3.5 

36:50:10.5 

39.51 ±0.03 

154 ±6 

6.33 

± 

0.07 

0.14 ±0.01 

18.0 ± 1.0 

27.0 ±6.0 

120.66 

6.78 

2 

NGG6207 

16:43:3.8 

36:50:3.4 

39.36 ± 0.03 

131 ±4 

6.15 

± 

0.06 

0.148 ±0.009 

13.9 ±0.9 

21.0 ±5.0 

0.51 

7.38 

1 

NGG1073 

2:43:38.3 

1:23:38.6 

39.22 ±0.03 

124 ±3 

6.04 

± 

0.05 

0.138 ±0.007 

13.4 ±0.5 

23.0 ±3.0 

1868.64 

6.39 

2 

NGG1073 

2:43:39.6 

1:23:43.3 

38.87 ±0.06 

87 ±3 

5.64 

± 

0.08 

0.16 ±0.01 

9.9 ±0.6 

23.0 ±5.0 

521.43 

6.59 

1 

NGG3423 

10:51:19.2 

5:51:8.5 

39.64 ± 0.02 

407 ±5 

7.03 

± 

0.03 

0.038 ±0.001 

14.5 ±0.3 

9.3 ±0.8 

17.46 

6.95 

2 

NGG3423 

10:51:18.9 

5:51:19.4 

39.37 ±0.02 

236 ±8 

6.54 

± 

0.06 

0.063 ± 0.004 

6.7 ±0.6 

< 4.0 ± 1.0 

16.14 

6.9 

1 

NGG428 

1:12:53.2 

0:56:47.5 

39.27 ±0.02 

167 ±2 

6.26 

± 

0.03 

0.093 ± 0.003 

10.8 ±0.4 

12.0 ± 1.0 

28.12 

6.83 

2 

NGG428 

1:12:53.2 

0:56:47.5 

39.03 ± 0.03 

155 ±4 

6.09 

± 

0.06 

0.08 ± 0.005 

5.7 ±0.6 

5.0 ± 1.0 

28.22 

6.83 

1 

NGG5334 

13:52:55.8 

-1:5:18.9 

38.93 ±0.01 

171 ±6 

6.11 

± 

0.06 

0.061 ± 0.004 

< 4.0 


57.57 

6.8 

2 

NGG5334 

13:52:52.3 

-1:6:46.8 

38.81 ±0.02 

172 ±7 

6.05 

± 

0.07 

0.053 ± 0.004 

< 4.0 


29.0 

6.83 

1 

NGG5112 

13:22:2.4 

38:43:15.2 

39.55 ± 0.04 

288 ±7 

6.75 

± 

0.06 

0.057 ±0.003 

11.6 ±0.5 

8.0 ± 1.0 

155.01 

6.77 

2 

NGG5112 

13:22:2.0 

38:43:15.0 

39.54 ± 0.04 

215 ±6 

6.56 

± 

0.06 

0.088 ± 0.005 

12.4 ±0.6 

11.0 ±2.0 

137.97 

6.78 

1 

NGG3403 

10:53:56.6 

73:41:27.7 

38.78 ± 0.05 

200 ± 10 

6.1 

± 

0.1 

0.041 ± 0.005 

9.2 ±0.8 

15.0 ±5.0 

13.86 

7.03 

2 

NGG3403 

10:53:58.5 

73:41:19.6 

38.71 ±0.05 

171 ±9 

6.0 

± 

0.1 

0.048 ± 0.005 

4.2 ±0.8 

< 4.0 ±2.0 

39.29 

6.82 

1 

NGG4639 

12:42:55.5 

13:14:53.2 

38.63 ± 0.08 

160 ±9 

5.9 

± 

0.1 

0.047 ±0.006 

10.8 ±0.8 

26.0 ±8.0 

50.94 

6.81 

2 

NGG4639 

12:42:49.1 

13:16:7.4 

38.33 ± 0.03 

238 ±6 

6.02 

± 

0.05 

0.0186 ±0.0009 

7.8 ±0.9 

16.0 ±5.0 

75.81 

6.8 

1 

NGG7241 

22:15:49.3 

19:13:49.7 

38.64 ± 0.03 

126 ±5 

5.76 

± 

0.07 

0.069 ± 0.005 

< 4.0 


130.24 

6.78 

2 

NGG7241 

22:15:49.7 

19:13:48.1 

38.54 ± 0.08 

106 ±5 

5.6 

± 

0.1 

0.08 ± 0.008 

8.0 ±2.0 

20.0 ± 10.0 

94.91 

6.79 

2 

NGG918 

2:25:49.6 

18:30:25.2 

38.58 ± 0.04 

87 ±3 

5.49 

± 

0.07 

0.111 ±0.007 

10.2 ±0.7 

34.0 ±8.0 

52.28 

6.81 

1 

NGG4324 

12:23:6.2 

5:15:1.7 

39.31 ±0.02 

165 ±5 

6.28 

± 

0.05 

0.1 ±0.005 

< 4.0 


5.74 

7.1 

2 

NGG4324 

12:23:5.1 

5:14:53.4 

37.97 ±0.06 

84 ±5 

5.2 

± 

0.1 

0.058 ±0.007 

< 4.0 


10.56 

7.06 

1 

NGG4389 

12:25:31.5 

45:43:5.6 

38.98 ± 0.02 

161 ±6 

6.09 

± 

0.07 

0.071 ± 0.005 

13.0 ± 1.0 

25.0 ±6.0 

12.86 

7.04 

2 

NGG4389 

12:25:31.4 

45:43:5.1 

38.59 ±0.01 

199 ±9 

6.03 

± 

0.07 

0.033 ± 0.002 

4.6 ±0.7 

< 4.0 ±2.0 

15.44 

6.94 

1 

NGG4498 

12:31:39.5 

16:51:22.2 

38.7 ±0.03 

120 ±3 

5.76 

± 

0.06 

0.079 ± 0.005 

13.3 ±0.7 

43.0 ±8.0 

30.39 

6.83 

2 

NGG4498 

12:31:41.4 

16:50:51.4 

38.6 ±0.06 

80 ±3 

5.44 

± 

0.09 

0.13 ±0.01 

5.7 ±0.6 

11.0 ±3.0 

70.1 

6.8 

1 

NGG2541 

8:14:40.4 

49:2:30.7 

39.48 ± 0.04 

160 ±3 

6.34 

± 

0.05 

0.127 ±0.006 

18.4 ±0.3 

29.0 ±3.0 

44.26 

6.81 

2 

NGG2541 

8:14:36.3 

49:3:7.6 

39.19 ±0.03 

119 ±3 

6.0 ±0.06 

0.143 ±0.009 

14.2 ±0.5 

28.0 ±5.0 

90.61 

6.79 

1 

NGG2500 

8:1:46.9 

50:44:3.3 

38.44 ± 0.04 

76 ±4 

5.3 ±0.1 

0.12 ±0.01 

7.7 ±0.6 

24.0 ±7.0 

214.74 

6.76 

2 

NGG2500 

8:1:59.7 

50:45:57.8 

38.4 ± 0.04 

99 ±2 

5.49 ±0.05 

0.075 ± 0.004 

9.4 ±0.8 

33.0 ±8.0 

20.71 

6.96 



